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New ZEUS data are added into the NLO QCD analysis of the ZEUS- JETS PDF fit. 
The addition of high-Q NC and CC e~p inclusive cross-section data improves the 
determination of the it-valence quark at high x. The addition of high-Q 2 CC e + p 
inclusive cross-section data improves the determination of the d- valence quark at high 
x. The addition of lower-Q 2 NC e + p inclusive cross-section data, run at three different 
proton beam energies, improves the determination of the sea and gluon PDFs at small 
x. The new PDF fit is called the ZEUS09 PDF fit. 



1 Introduction 

The kinematics of deep inelastic lepton hadron scattering (DIS) is described in terms of 
the variables Q 2 , the negative invariant mass of the exchanged vector boson, Bjorken x, 
the fraction of the momentum of the incoming nucleon taken by the struck quark (in the 
quark-parton model), and y which measures the energy transfer between the lepton and 
hadron systems. The double differential cross-sections for the neutral current (NC) process 
with lepton polarization P are given by, 

d 2 a{e ± p) _ W r ■ ±1 ± ,p o,P m 

dxdQ 2 -^ [Ho+PHpl > H ^- Y+F2 ~ yF t TY - xF z W 

where, Y± = 1 ± (1 — y) 2 , and, at LO in QCD, the structure functions F° ' P and xF®' P are 
directly related to quark distributions by 

F° 2 P = £ *(* + - qi )A«> P , xF° 3 P - E - «)^ P ( 2 ) 
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The coefficients A®, for unpolarised beams are given by 

A^Q 2 ) = e 2 - 2e l v l v e Pz + (v 2 e + a 2 e )(v 2 + af)P 2 (3) 

B°(Q 2 ) = ~2e iai a e P z + AaiViV e a e Pi (4) 
The coefficients for the polarisation terms are given by 

A p = 2e l a e v l P z - 2a e v e {v 2 + a 2 )P%, (5) 
Bf = 2e i a i v a P z -2a i Vi(v* + al)P%. (6) 
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The term in Pz arises from 7Z interference and the term in P| arises purely from Z° 
exchange, where Pz accounts for the effect of the propagator relative to that of the 
virtual photon, and is given by 

P Q2 1 (7) 
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The other factors in the expression for A and B are the charge, e,, NC electroweak vector, Vi, 
and axial- vector, cij, couplings of quark i and the corresponding NC electroweak couplings 
of the electron, r, ,a e . Equations I5|6I show that polarization effects are only important at 
high Q 2 . 

Unpolarized HERA data have been used in fits to determine Parton Distribution Func- 
tions (PDFs). For low x, x < 10 -2 , F® is sea quark dominated and its Q 2 evolution, as 
predicted by QCD, is controlled by the gluon contribution, such that HERA data provide 
crucial information on low- a; sea-quark and gluon distributions. At high Q 2 , the structure 
function xF® becomes increasingly important, and gives information on valence quark dis- 
tributions. The charged current (CC) interactions also enable separation of the flavour of 
the valence distributions at high-x, since their (LO) cross-sections are given by, 
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Parton Density Function (PDF) determinations are usually global fits [HOE], which use 
fixed target DIS data as well as HERA data. In such analyses, the high statistics HERA 
NC e + p data have determined the low-x sea and gluon distributions, whereas the fixed 
target data have determined the valence distributions. Now that high-Q 2 HERA data on 
NC and CC e + p and e~p inclusive double differential cross-sections are available, PDF fits 
can be made to HERA data alone, since the HERA high Q 2 cross-section data can be used 
to determine the valence distributions. This has the advantage that it eliminates the need 
for heavy target corrections, which must be applied to the v-Fe and fiD fixed target data. 
Furthermore there is no need to assume isospin symmetry, i.e. that d in the proton is the 
same as u in the neutron, since the d distribution can be obtained directly from CC e + p 
data. 

The ZEUS-JETS PDF fit [5] was an NLO QCD fit in the DGLAP formalism to ZEUS in- 
clusive cross-section data and jet production data from HERAT. The PDFs were parametrized 
at Qq = 7GcV 2 by the form xf(x) = pix P2 (l — x) P3 (l +P4X), using 11 free parameters. 
Predictions for the cross-sections were made by evolving the PDFs to the Q 2 values of 
the measurements and convoluting them with coefficient functions, calculated in the general 
mass variable flavour-number scheme of Thorne and Roberts |13j , to produce structure func- 
tion predictions. Predictions for jet cross-sections were made using NLO programmes [Sl[7]- 
In evaluating the uncertainty on the PDF parameters and the cross-section predictions which 
derive from them full acount is taken of correlated experimental uncertainties using the Off- 
set method. However the determinations of the valence PDFs from HERA-I data are not as 
accurate as those from global fits because of poor statistics at high-x. The addition of data 
from HERA-II changes this situation. In order to assess the impact of this new data the 
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Figure 1: ZEUS NC (left-hand side) and CC (right-hand side) e~p data from HERA-II 
running with polarised beams. The predictions of the ZEUS09 fit are superimposed 

fit formalism has not been changed. This paper describes the improvement to the valence 
PDFs from adding ZEUS e~p NC [Sj from the 2005-6 running period, e~p CC j9] data from 
the 2004-6 running period, and e + p CC data [10] from the 2006-7 running period, into the 
ZEUS-JETS fit. 

The polarization of the HERA-II data can also be exploited to measure electroweak 
couplings. The CC cross-sections give information on the propagator mass and the weak 
coupling. The NC cross-sections give information on the quark couplings to Z°. A prelim- 
inary model independent extraction of these parameters was given in Ref [11] and will be 
updated when the full HERA-II data set is available. In the present paper we compare our 
data to the electro- weak predictions of the Standard Model. 

Most of the data collected at HERA were for collisions between e of 27.5GcV with 
protons of 920GeV beam energies. During the final running period, HERA provided lower 
Q 2 data (24 < Q 2 < 110 GeV 2 ), with modified trigger conditions, collected at three different 
proton beam energies (920, 575, 460GeV). These data access high-j/ and have been used to 
measure the longitudinal structure function Fl P~2]. At low- a;, NLO QCD in the DGLAP 
formalism predicts that this structure function is strongly related to the gluon PDF. The 
reduced cross-section data from these runs have also been added into the ZEUS- JETS PDF 
fit and this provides an improved determination of the sea and gluon PDFs at small-x. The 
new fit including all these new data is called the ZEUS09 PDF fit. 

2 Results 

Fig Q] (left-hand side) shows ZEUS data on NC e~p double differential cross-sections from 
the 2005-2006 running period. There are 99p6 _1 of negatively polarised data (P e = —0.27) 
and 71pb~ 1 of positively polarised data (P e = +0.29). This figure (right-hand side) also 
shows ZEUS data on CC e~p double differential cross-sections from the 2004-2006 running 
period. There are 104pfo~ 1 of negatively polarised data (P e = —0.27) and 71pfr~ J of positively 
polarised data (P e = +0.30). Fig [2] (left-hand side) shows ZEUS data on CC e + p double 
differential cross-sections from the 2006-2007 running period. There are bGpb^ 1 of negatively 
polarised data (P e = —0.36) and 76p& -1 of positively polarised data (P e = +0.33). The 
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Figure 2: ZEUS CC e + p data from HERA-II running with polarised beams (left-hand side). 
ZEUS NC e + p data from 2007 HERA-II running with three different proton beam energies 
(right-hand side) . The predictions of the ZEUS09 fit are superimposed 

same figure (right-hand side) shows ZEUS data on NC e + p double differential cross-sections 
for three proton beam energies corresponding to, s/s = 318,251,225 GeV, from the 2007 
running period. There are 44.5, 7.1, 14. Opb^ 1 of data for each beam energy respectively. 
After the addition of these new data the 577 data points in the ZEUS- JETS PDF fit has 
increased to 1060 data points in the new ZEUS09 PDF fit. The results of this ZEUS09 fit 
are superimposed on the data in all these figures. The \ 2 P er degree of freedom of this fit is 
0.97. The agreement of all of the polarised data with the fit is a confirmation not only of the 
validity of the NLO QCD in the DGLAP formalism but also of the electroweak predictions 
of the Standard Model in a space-like process. 

Fig. [3] (left-hand-side) illustrates the further impact of adding these new data by com- 
paring the fractional uncertainties of the PDFs extracted from a fit including these data 
with those extracted from the ZEUS- JETS PDF. The improvement in the u- valence quark 
at high x comes from the addition of the NC and CC e~p data from the 2004-2006 running. 
Both these processes are u quark dominated at large x. The d-valence quark uncertainty 
is also reduced significantly at large x. This improvement derives from the CC e + p data 
which are d quark dominated at large x. Finally the improvement in the low- a; gluon and 
sea PDFs come from the addition of the NC e + p lower Q 2 data run at three different beam 
energies. This final addition of data completes the ZEUS09 PDF fit. The \ow-x scale of 
the figure is extended to illustrate the improvement which is expected from adding these 
high-y data. The PDFs extracted from the ZEUS09 fit are also compared to those of the 
ZEUS- JETS fit in Fig [3] (right-hand side). The central values of the fit are very compatible 
with the ZEUS- JETS fit, but the gluon PDF is a little steeper indicating the impact of the 
low-Q 2 data on the gluon PDF. 

3 Summary 

The inclusion of high-Q 2 NC and CC e~p data into the ZEUS- JETS PDF fit results in an 
improved determination of the it-valence PDF. The further inclusion of high-Q 2 CC e + p 
data results in an improved determination of the li-valencc PDF. These data were run with 
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Figure 3: Left hand side: The fractional uncertainties of the ZEUS09 PDFs compared to 
those of the ZEUS- JETS PDFs. Right hand side: The PDFs extracted from the ZEUS09 
fit compared to the ZEUS- JETS fit at Q 2 = lOGcV 2 . 

polarised lepton beams and the cross-section data for the different polarisations provides 
a spectacular confirmation of the Standard Model electroweak predictions in a space-like 
process. Finally the inclusion of lower-Q 2 data run at three different proton beam energies 
yields an improved determination of the low- a; sea and gluon PDFs. The new fit including 
all these data is called the ZEUS09 PDF fit pp. 
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